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PREFACE 

The  development  of  automatic  radar  detection  systems  has 
focused  attention  on  the  effects  on  the  performance  of  a  radar 
receiver  of  the  interfering  pxilses  from  other  radars.  It  is 
essential  to  the  successful  performance  of  aan  automatic  system 
that  the  false  indications  caused  by  interference  be  almost 
completely  suppressed.  This  paper  is  a  mathematical  investiga- 
tion of  one  type  of  radar  receiver  which  discriminates  against 
false  target  indications.  A  double  threshold  radar  receiver  is 
found  to  be  capable  of  controlling  the  false  target  interference. 
The  double  threshold  receiver  is  shovm  to  have  a  slightly  lower 
radar  range  capability  than  a  normal  receiver. 

Sincere  thanKs  are  due  Mr.  John  Mallet t  and  Dr.  i'eter  bwerling 
of  the  RAND  Corporation,  Santa  Monica,  California,  for  their 
guidance  and  assistance  in  the  pursuit  of  this  investigation. 
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CHAPTER  I 
INTRODUCTION 

The  continuing  development  of  radar  systems  in  recent  years  has  led 
to  an  increasing  usefullness  of  these  systems  for  search  and  for  control 
functions.  The  operational  utilization  of  radar  by  the  armed  forces, 
even  though  high  at  the  end  of  world  War  II,  has  continued  to  grow  ever 
since.  Operational  planners  have  called  for  greater  and  greater  numbers 
of  radars  to  meet  the  service  requirements,  ^t  with  this  growth  in 
numbers  of  radars  the  effect  that  they  would  have  upon  each  other  has 
been  rather  generally  ignored.  A  basic  analysis  of  the  numbers  of  radars 
involved  in  current  plans  indicates  that  a  very  high  signal  density  is 
developing  in  all  operating  areas.  A  radar  receiver  in  this  operating 
situation  will  intercept  the  pulses  transmitted  by  many  other  radars  in 
the  area  and  will  have  to  contend  with  the  large  numbers  of  false  indica- 
tions which  these  interfering  pulses  will  cause. 

The  presence  of  randomly  occurring  interfering  pulses  intercepted 
by  a  manned  search  radar  are  presently  not  a  primary  limitation  on  the 
performance  of  the  radar  system.  These  interfering  pulses,  tne  trans- 
mitted signals  of  other  radars,  appear  on  the  PPI  scope  of  a  search  radar 
as  false  targets j  however,  in  nearly  alx  cases  proper  training  of  the 
radar  operator  permits  him  to  evaluate  them  correctly,  and  the  overall 
system  (radar-plus -operator)  performs  properly.   In  some  dense  signal 
environments  the  training  necessary  to  permit  the  operator  to  identify 
false  targets  caused  by  interfering  signals  ma^/   become  lengthy.  However, 
he  eventually  is  able  to  carry  out  his  role  as  the  essential  evaluation 


oleaent  in  the  radar  system.  In  some  locations  in  England  today,  vrtiere 
the  density  of  radars  is  relatively  high,  the  numbers  of  interfering 
pulses  are  so  great  that^  a  new  operator  in  a  GCI  radar  installation  re- 
quires approximately  one  week  of  training  to  become  familiar  with  the 
many  false  targets  presented  by  the  transmissions  of  other  radars. 

Ifith   the  advent  of  more  nearly  automatic  systems,  the  presence  of 
false  targets  beccanes  a  primary  limitation  on  system  performance.  With 
new  air  search  radar  systems  and  with  guided  missile  systems  the  trend 
is  to  remove  the  operator  from  the  system  and  to  substitute  automatic 
computers  of  various  degrees  of  complexity.   In  these  systems  the  random- 
ness in  space  and  in  time  of  the  false  targets  makes  it  extremely  diff- 
icult if  not  impossible  to  handle  them  properly  in  the  computers,  inter- 
fering pxilses  can  cause  the  performance  of  these  automatic  and  semi- 
automatic systems  to  break  down  completely.  The  computers  cannot  eval- 
uate large  numbers  of  false  targets  correctly,  and  the  system  approaches 
such  a  state  of  saturation  that  an  actual  target  would  very  probably  be 
missed.  The  large  numbers  of  radars  now  coming  into  xise  in  oxir  armed 
services  are  causing  an  ever-increasing  density  of  signals  which  the 
automatic  systems  must  contend  with.  Consideration  of  present  radars 
indicates  that  a  search  radar  in  the  x-band  would  now  be  confronted  with 
an  interference  density  level  of  several  thousand  pulses  per  second. 
Acquisition  of  new  radar  systems  will  steadily  increase  this  figure. 

The  development  of  a  radar  receiver  that  can  eliminate  large  numbers 
of  interfering  pulses  without  a  significant  loss  in  its  actual  target 
detection  capability  becomes  a  matter  of  paramount  concern,   if  this  can 
be  accomplished,  then  the  data  actually  fed  to  the  computers  can  be 
handled  correctly  resulting  in  an  overall  efficient  and  accurate  system 


performance.  This  report  is  an  investigation,  by  use  of  a  mathematical 
model,  of  a  method  for  handling  interfering  pulses  in  the  radar  receiver. 

One  characteristic  of  the  interfering  pulses  vrtiich  can  be  useful  in 
distinguishing  them  from  target  echoes  is  their  randomness  of  occurrence 
in  space  and  time.  As  the  rotating  antenna  of  a  search  radar  scans  past 
a  target  all  of  the  target  echoes  occur  at  very  nearly  the  same  range. 
With  most  present  air  search  radars  the  antenna  rotation  rate,  the  beam 
width,  and  the  pulse  repetition  frequency  are  such  that  from  ten  to  fifty 
echo  pulses  can  be  expected  to  be  returned  f^om  a  target  in  one  scan  of 
the  antenna.  During  this  short  interval  of  time  (approximately  0.05  to 
0.1  seconds)  the  target  range;  cannot  change  appreciably  (in  0.1  second 
a  1200  knot  aircraft  travels  approximately  70  yards).  The  interfering 
pulses  which  may  be  intercepted  by  the  radar  in  this  same  time  interval 
occur  at  random,  and  the  probability  that  more  than  a  very  few  occur  at 
times  corresponding  to  one  range  in  the  search  radar  becomes  very  small. 
It  is  proposed  that  this  difference  between  interfering  pulses  and  actual 
t^arget  echoes  be  utilized  in  a  radar  receiver  system  designed  to  suppress 
the  false  targets. 

A  receiver  with  this  capability  would  conceivably  be  designed  with 
its  effective  search  range  divided  into  discrete  range  intervals.   If 
each  of  these  range  intervals  corresponded  to  the  pulse  duration,  then 
the  returned  target  echoes  >.ould  occur  in  one  interval  or  in  two  adjacent 
intervals.  A  target  would  be  announced  only  when  the  number  of  signals 
occurring  in  one  range  interval  in  the  time  equivalent  to  the  antenna 
scanning  past  a  target  exceeded  some  pre assigned  number.  This  number 
would  be  established  high  enough  to  miniiiiize  the  probability  of  a  false 
target  indication  from  randomly  occurring  interfering  pulses  while  still 
permitting  an  actual  target  to  be  recognized. 


CHAPTER  II 
SIGNAL  DETECTION  THEORIES 

Before  proceeding  with  a  detailed  analysis  of  the  performance  of  a 
radar  receiver  of  this  type  it  would  be  well  to  review  the  work  on  radar 
signal  detection  theory  that  has  been  accomplished  in  recent  years, 
workers  at  the  Bell  Laboratories  and  at  the  RAND  Corporation  have  identi- 
fied the  statistical  nature  of  the  noise  in  a  radar  receiver  and  have 
employed  the  laws  of  statistics  to  develop  a  theory  for  the  probability 
of  detection  of  a  target. 

The  basis  of  all  of  the  theories  on  the  statistical  nature  of  signal 
detection  is  the  work  of  s.  0.  Rice/  5  /•  He  has  defined  the  probability 
distribution  fxmction  for  noise  and  has  developed  the  distribution  func- 
tions that  deecribe  t,he  effect  of  filtering. 

Both  the  thermal  noise  voltage  across  a  resistor  and  the  noise  volt- 
age due  to  the  shot  effect  in  a  vacuum  tube  approach  a  normal  distri- 
bution vrtien  the  number  of  electrons  involved  per  second  in  the  processes 
tends  toward  infinity.   In  practice,  it  may  usually  be  assumed  that  the 
total  noise  voltage  between  any  two  points  due  to  any  combination  of 
thermal,  shot,  and  cosmic  noise  sources  can  be  represented  hy   the  distri- 
bution function 
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P  =  "pTtW^   e  (1) 

where  ^fi  is  the  mean  square  value  of  the  noise  voltage.  This  distri- 
bution is  valid  provided  all  elements  involved  in  the  composition  of  the 
total  noise  voltage  have  been  linear. 

If  this  noise  is  now  passed  through  a  linear  filter  of  center  fre- 
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quency  f^   ,  having  a  pass  band  which  is  narrow  compared  to  f^  ,  the 
output  will  have  an  envelope  which  has  a  probability  density  function 

R    -  R^/2  'f  » 
Pr-^e  (2) 

where  R  is  the  amplitude  of  the  envelope.  (Pc  is  the  mean  square 

noise  voltage,  given  by 

^(f)       df  (3) 

w(f )  is  the  power  spectrum  of  the  filter. 

If  the  input  to  a  filter  consists  of  a  sine  wave  of  frequency  fj^, 
as  well  as  noise,  then  the  probability  density  function  of  the  output 
envelope  is  o       o 

(rI^p^)         /     \ 

p  —     r      --  .1     -      •  /    pp      \ 

PzTjT    e         ^f^         lo[^j  R>  0  (U) 

Pro  R<i^  0 

where  p  is  the  amplitude  that  the  sine  wave  would  have  at  the  output 
of  the  filter  in  the  absence  of  noise  and  lo  is  the  modified  Bessel 
function  of  the  first  kind. 

The  envelope  of  the  output  has  a  correlation  time  which  is  approxi- 
mately equal  to  the  reciprocal  of  the  bandwidth  of  the  filter.  Thus, 
it  is  improbable  that  the  envelope  vrill  change  by  an  appreciable  per- 
centage in  times  much  less  than  the  correlation  time,  but  it  is  quite 
probable  that  it  will  change  by  a  goodly  percentage  in  times  large 
compared  with  the  correlation  time.  Jt  is  therefore  assumed  that  values 
of  the  envelope  1/ A^  seconds  apart  are  independent,  where  ^f  is  the 
bandwidth  of  the  filter. 

Starting  with  these  basic  expressions,  J.  J.  Marcum  of  the  RAND 
Corporation  /  2  /  and  R.  L.  Kaplan  of  Bell  Laboratories  /  1  /  have 


independently  developed  theories  of  signal  detection  by  pulsed  radar. 
These  two  have  used  slightly  different  matheaiaticai  definitions  for 
the  process  of  detection.  For  this  reason,  the  mathematical  ateps 
leading  to  the  final  result  are  different.  The  final  results  in  the 
two  developments  can  be  shown  to  be  equivalent. 

Marcura  assumes  that  the  radar  receiver  is  a  linear  integrator. 
In  this  concept  the  receiver  linearly  adds  the  voltage  output  of 
N  samples  from  the  rectii'ier.  Because  of  the  time  elapsing  between 
sampling,  these  samples  are  considered  to  be  independent.  The  samples 
are  considered  to  be  either  of  signal-plus-noise  or  of  noise  alone.  If 
the  sum  of  the  N  samples  exceeds  the  bias  level  calculated  from  the 
probability  density  function  for  N  variates  of  noise  alone,  then  a 
signal  is  said  to  be  present. 

The  integrator  may  take  the  sum  of  the  squares  of  the  N  variates, 
or,  in  general,  the  sum  of  N  variates  where  each  variate  has  been 
processed  by  some  general  function,  as  long  as  the  same  weight  is 
applied  to  each  variate,  the  integrator  is  called  linear. 

The  best  possible  function  to  produce  integrable  equations  in 
the  development  of  the  detection  theory  is  that  given  by  a  square  law 
rectii'ier.  Both  M^rcun  and  Kaplan  show  that  linear  rectification 
does  not  give  very  different  results  (  about  0.2  db). 

In  discussing  Marcum's  results  the  following  terms  will  be  used: 
Pdjj    ~  probability  of  detectionj  the  probability  that 
the  sum  of  n  samples  of  signal-plus-noise  will 
exceed  the  bias  level. 
Y     -  integrator  output  for  the  sum  of  N  variates. 
Yb     -  bias  level. 
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X     -  signal-to-noise  power  ratio  of  the  input  signal, 
ll     -  probability  of  a  false  target;  the  probability 
that  the  sum  of  N  samples  of  noise  alone  will 
exceed  the  bias  level. 

Since  the  receiver  noise  is  normally  distributed,  there  is  a  finite 
probability  that  the  sum  of  N  noise  variates  alone  will  exceed  the  bias 
level,  Yb>  J^o  matter  how  high  this  may  be  set.  This  will  result  in  a 
false  target  indication,  or  a  false  alarm.  In  Marcum's  mathematical 
model  of  a  radar  receiver  this  bias  level  is  established  at  a  level 
which  will  give  a  desired  probability  of  false  alarm  from  noise  alone. 

This  false  alarm  probability  can  be  related  to  the  time  interval 
between  false  alarms.  As  explained  above,  values  of  the  receiver  out- 
put envelope  l/^  f  seconds  apart  are  independent.  Since  the  receiver 
^©ise  is  present  at  all  times,  in  a  time  interval  t  there  are  T*^  f 
independent  chances  for  a  false  alarm  to  occur  from  noise  alone.  Thus 
the  probability  of  the  occurrence  of  at  least  one  false  alarm  is  a 
function  of  time.  Marcum  has  chosen  to  define  the  false  alarm  time 
as  the  time  in  which  the  probability  is  one-half  that  noise  alone  will 
not  exceed  the  bias  level.  He  has  defined  the  false  alarm  number  as  the 
T»^  f  independent  chances  for  a  false  alarm  in  this  false  alarm  time,  T. 

With  this  description  of  the  mathematical  model,  Marcum's  results 
can  be  summarized: 
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T,     .    ^    in  2  4  (0.693)  M 

where  n  is  the  false  alarm  niimber  defined  above     (n  =  T»A   f). 


where   ll^i   ,  N  -  ij  is  the  incoaplete  gamma  function.  Tables  of 
I  (u,  p)  have  been  compiled  by  Pearson  /  \xj     * 

In  practice  n  and  N  are  given  and  the  desired  result  is  the  bias 
level  Yb-  Since  the  equation  cannot  be  solved  directly  f or"  y^,  com- 
putations are  facilitated  if  curves  are  first  plotted  showing  the 
general  relationship  of  equations  (i?)  and  (6). 

Knowing  the  required  bias  level  for  a  given  false  alarm  interval, 
the  cumulative  probability  of  detection  for  n  variates  of  signal -plus- 
noise  is 
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where  I„ _•,  (2  ^NTi)  is  the  modified  Bessel  function. 

This  integral  is  not  solvable  directly  in  terms  of  well-known 
functions.  For  N  =  1  (the  probability  that  one  signal -pi us -noise 
pulse,  with  signal-to-noise  ratio  x,  exceeds  the  bias  level  Y^)  the 
equation  has  been  solved  by  machine  computation  at  the  RAND  Corporation. 
Since  this  function  will  be  Uiied  in  later  developiaents,  the  solution 
is  presented  graphically  in  Figure  i?U  for  various  values  of  Y^* 

For  N";^!  the  solution  to  the  detection  integrad  has  been  shown 
to  be 

Pd^  =  ^  "  5y^  (2N-1,  N-1,  /n?)  (b) 

where  Tg  (m,n,r)  =  2  r  ""^*^  e'^^C^  t'^""  ^-t^  i^(2rt)  dt 

is  the  incomplete  Toronto  function.  Curves  of  this  function  are  avail- 
able in  reference  (2 ) . 


Peter  swerling,  also  of  the  RAND  Corporation,  has  extended 
Marcum's  basic  theory  in  two  reports,  in  the  first  /  6  /  he  has 
considered  the  effects  of  fluctuation  in  the  amplitude  of  target 
echo  signals  on  the  probability  of  target  detection,   in  general, 
two  situations  were  investigated:   (1)  the  target  echo  signals  would 
have  constant  amplitude  during  the  time  for  the  radar  antenna  to 
scan  past  a  target  but  would  fluctuate  independently  from  scan  to 
scan;   and  (2)  the  amplitude  of  the  target  echo  signals  would  fluctu- 
ate independently  from  pulse  to  p\ilse. 

For  scan  to  scan  fluctuations  sterling  found 

Y>,_ 

1  -  1  .       Pdi  =  6   '  *  "  (9) 
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where  I  (u,  p)  is  the  incomplete  gamma  function  and  x  is  the  mean  over 
all  values.  The  target  range  at  vriiich  there  was  a  90<?  probability 
of  detection  with  scan  to  scan  target  echo  fluctuations  was  found  to 
be  about  3/5  of  that  found  for  constant  target  echoes.  The  range  for 
50%   probability  of  detection  with  scan  to  scan  fluctuations  was  about 
9/10  of  the  range  for  constant  echoes. 

For  pulse  to  pulse  fluctuations  sterlings  results  were 


The  range  for  90^  probability  of  detection  with  pulse  to  pulse 
target  echo  fluctuations  was  about  19/20  of  that  for  constant  target 
echoes,  and  the  range  for  50$  probability  of  detection  was  approximately 
equal  to  the  range  for  constant  echoes. 

In  sterling's  second  report  /  7  /  a  mathematical  analysis  was  made 
of  a  so-called  "double  threshold"  receiver,  considering  only  interfer- 
ence from  receiver  noise.  This  is  the  type  of  receiver  described 
earlier  as  having  possible  advantages  in  the  suppression  of  false 
alarms  from  interfering  pulses.  The  first  threshold  is  the  bias 
level  Yb»  The  receiver  counts  the  number  of  pulses  which  exceed  this 
bias  level.  A  target  is  declared  to  be  present  if  this  number  is  at 
least  equal  to  some  fixed  number  M.   If  there  are  N  target  echoes 
per  scan,  then  M^' N.  The  number  M  is  called  the  second  threshold. 

If  p^  represents  the  probability  that  a  single  signal-plus- 
noise  pulse  exceeds  the  bias  level  y^  (as  defined  in  equation  (7)  and 
illustrated  graphically  in  figure  Sh)   then  swerling's  results  can  be 
stated 

_.N-,.  N         N-i     i 
^N-f-'^i   (^-Px)     Pi  (12) 

where  c^  stands  for  the  binomial  coefficients. 

Swerling  found  that  this  method  has  a  0.5  to  1.0  db  loss  in 
sensitivity  as  compared  to  a  normal  integrating  receiver.  He  also 
found  that  for  any  n  there  is  an  optimum  value  for  m/n  which,  in 
the  range  of  interest  to  present  search  radar  systems,  varied  from 
m/n  z  0.6   for   N  -  10   to  •  m/n  r  O.U   for    n  z  UO. 

E.  L.  Kaplan  has  assumed  the  following  detection  procedure.  The 
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rectified  output  is  sampled  at  n  different  instants  of  time  which 
are  far  enough  apart  so  that  the  values  of  the  noise  are  independent. 
If  the  average  output  at  these  instants  exceeds  the  average  noise 
level  Y  by  an  amount  ky  or  more,  a  signal  is  announced.  The  noise 
envelope  at  the  output  of  a  narrow  band  filter  has  a  Rayleigh  distri- 
bution (equation  (2)  ),  The  average  of  n  such  values  (assuming  in- 
dependence) has  a  chi-square  (^^)  distribution  with  2n  degrees 
of  freedom.  This  distribution  is  manipulated  to  obtain  the  desired 
expressions  for  the  res\ilts. 

m  a  discussion  of  these  results,  u  is  the  standard  normal 
deviate.  And,  using  the  previous  notation, 


r„= 


=    Pr(u^p)  =  \         _e jiu 


vAp    JTrf 


Pd  =  Pr(u>up)  =  j         ^ _ .         du 

For  the  steady  target  echo  case,  Kaplan's  results  give  the  signal-to- 
noise  ratio  required  to  achieve  the  desired  probability  of  detection 
with  the  desired  false  alarm  probability. 


^ 


f   --'■■■■   --      (X3) 


N 


Approximations  in  the  deprivations  introduce  some  inaccuracy  at  n  -  1. 
curves  of  x  as  a  function  of  n  are  plotted  in  reference  (1)  for 
the  steady  target  echo  case  and  for  the  fluctuating  target  echo  cases. 
KaTDlan  discusses  the  double  threshold  receiver  and  concludes  that  it 
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requires  a  signal-to-noise  ratio  about  l.U  db  above  that  Triiich  sufficed 
for  equal  performance  using  averages. 

All  of  this  work  on  signal  detection  theory  has  been  an  application 
of  the  stat^tical  theory  of  receiver  noise  to  the  pulsed  radar 
situation.  The  probability  of  a  false  target  indication  and  the 
false  alarm  number  have  been  based  on  noise  considerations  alone. 
The  effects  of  interfering  pulses  from  other  radars  have  not  been 
included,  although  they  will  obviously  increase  the  numbers  of  false 
alarms  appreciably. 

A  very  general  look  at  the  semi-automatic  radar  systems  coming 
into  operation  today  indicates  a  desired  false  alarm  time  (false 
alarms  from  noise  and  from  interfering  pulses  combined)  in  the  order 
of  a  few  minutes.  False  alarms  occurring  at  intervals  of  one  minute 
or  greater  can  be  handled  satisfactorily  in  the  automatic  computers. 
The  double  threshold  receiver  will  be  analyzed  'for  false  alarm 
times  of  one  minute  and  of  ten  minutes. 

preliminary  analyses  of  the  numbers  of  interfering  pulses 
occurring  in  dense  signal  enviornments  indicate  that  a  radar  in 
any  one  of  the  present  operating  bands  must  be  equipped  to  handle 
interference  amounting  to  thousands  of  pulses  per  second.  The 
double  threshold  receiver  will  be  analyzed  for  interference  levels 
of  10,000,  30,000,  and  i>0,000  pulses  per  second. 
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CHAPTER  III 
THE  DOUBLE  THRESHOLD  RECEIVER  -  NON-FLUCTUATING  TARGET  ECHOES 

A  double  threshold  receiver  counts  the  number  of  pulses'  whose 
amplitude  exceeds  a  preassigned  bias  level  Yfe  (first  threshold). 
Beyond  this  first  threshold  the  pulse  amplitude  is  not  considered. 
If  the  number  of  pulses  counted  in  one  range  gate  in  N  successive 
radar  repetition  periods  exceeds  a  preassigned  number  M  (second 
threshold),  a  target  is  declared  to  be  present. 

A  false  alarm  occurs  when  either  a  burst  of  noise  or  an  inter- 
fering pulse  occurs  in  one  range  gate  on  m  of  N  successive  repetition 
periods.  To  calculate  the  probability  of  a  false  alarm  it  is  necessary 
first  to  determine  the  probability  of  the  single  occurrence  of  a 
noise  pulse  or  an  interfering  pulse.  This  sixj^gle  occurrence  probability 
is  then  used  in  a  partial  sum  of  the  binomial  expansion  from  u     to 
N  to  arrive  at  false  alarm  probability. 

in  this  analysis  the  envelope  of  the  output  is  assumed  to 
have  a  correlation  time  which  is  approximately  equal  to  the  range 
gate  width.  Thus  a  noise  burst  occurring  in  a  range  gate  can  be 
considered  as  occupying  the  entire  range  gate,  and  values  of  the 
envelope  in  adjacent  range  gates  can  be  considered  to  be  independent 
as  far  as  random  noise  is  concerned. 

All  pulses,  both  target  signals  and  interfering  pulses,  are 
assumed  to  be  of  duration  equal  to  the  range  gate.  An  interfering 
pulse  occurring  at  randcan  in  the  repetition  period  will  not  coincide 
exactly  with  an  established  range  gate  but  will  overlap  so  that  a 
portion  of  the  piilse  is  in  each  of  two  successive  range  gates. 
Therefore,  each  pulse  is  considered  to  "ring"  two  adjacent  range  gates. 
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As  a  preliminary  step  in  investigating  false  alarms  from 
interfering  pulses  and  receiver  noise,  it  is  convenient  to  calculate 
the  effect  that  this  double  threshold  receiver  will  have  on  inter- 
fering pulses  alone*  The  following  additional  terms  will  be  used; 
p,-j  X  -  the  probability  that  one  or  more  pulses  will 

occur  in  one  range  gate. 
pL^    -  the  probability  that  pulses  will  occur  in  one 
range  gate  on  M  or   more  of  N  successive 
repetition  periods  (a  false  alarm). 
(J^       -   number  of  discrete  range  gates  per  second. 
T)  -  range  gate  width,  s  l/^  • 

S     -  average  number  of  random  interfering  pulses 
per  second. 
Assuming  a  random  arrival  of  interfering  pulses  (a  poisson  Model), 
the  probability  that  one  or  more  piiLses  will  occur  in  a  chosen  range 
gate  is 

P(l)  =  1  -  e    '-'         =  1  -  e    C  (lu) 

Since  in  cases  considered  here  2riS«!l  >  the  exponential  can  be 
represented  by  the  first  two  terms  of  its  series  expansion,  and 
the  desired  probability  becomes 

P(l)  =  \S  (15) 

Twice  the  range  gate  width  is  used  in  this  expression  since  a  pulse 
is  assumed  to  ring  two  adjacent  range  gates. 

The  probability  of  a  false  alarm  in  the  chosen  range  gate  becomes 

N   M         ^  N  -  i 

Pm  =  Z:  q      P(l)   (1  -  P(l)  )  (16) 


lU 


N 


V~   N       i  N-i 

Pm  =2-  C  (2-y^S)    (1  -  27|S)  (17) 

Figures  33  through  Ul  are  working  curves  that  show  the  partial  sum 
of  the  binomial  expansion  J^  7   from  M  to  N  for  N  =  10,  N  «  20, 
and  N  =  30.  After  p., v  Uas  been  found  from  the  radar  characteristics 
and  from  the  density  of  interference,  the  probability  of  false  alarm  in 
a  chosen  range  gate,  Py,  can  be  calculated  (equation  16  or  17)  by  use 
of  these  figures. 

The  receiver  is  assumed  to  have  <y        discrete  range  gates  per 
second,  or  y^]^ discrete  range  gates  in  each  repetition  period. 
A  double  threshold  radar  receiver  has  /f    possibilities  for  a  false 
alarm  in  a  time  interval  equivalent  to  N  repetition  periods. 
The  probability  of  false  alarm  in  this  interval  is 


^fa 


=  1  -  e"(  /4n)'^  (18) 


Figure  I  shows  the  effect  of  the  second  threshold  on  the  probability 
of  false  alarm,  from  interfering  pulses  for  N  =  10,  N  z   20,  and  n  =  30. 
In  this  figure  the  number  of  interfering  pulses  per  second  (s) 
multiplied  by  the  gate  width  (  ^y^  )  is  equal  to  0.05  {Sy)    -  O    -   'OS). 
This  corresponds  to  i>0,000  interfering  pulses  per  second  for  a  radar 
with  a  one  microsecond  gate  width. 

Although  this  figure  demonstrates  the  false  alarm  situation  over 
a  very  short  period  of  time,  it  is  useful  here  to  show  the  effect  of 
veu^ring  the  second  threshold,  M>  for  each  value  of  N.  Por  each  N  there 
is  a  very  definite  lower  value  for  jj  below  which  the  double  threshold 
technique  offers  no  advantages  in  eliminating  false  alarms.  For  values 
of  M  above  this  lower  limit  the  probability  of  false  alarm  decreases 
very  rapidly  with  increasing  m. 
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It  is  now  necessary  to  relate  the  probability  of  false  alarm  to 
a  longer  interval  of  time.   In  this  connection,  use  will  be  made  of  the 
false  alarm  time,  Tfa  »  previously  defined  by  Marcum  (i.e.,  the  time 
in  which  the  probability  is  one-half  that  a  false  alarm  will  not  occur). 

The  probability  of  at  least  one  false  alarm  in  time  T  is 


Pfa(T)  =  1  - 


^l-%l)  •  (^-%2)  '  (}-\3^ (^-V^ 


(19) 


where  pw;j_  is  the  probability  that  a  false  alarm  will  occur  in  the  first 
range  gate,  etc.;  and  r  ,  the  last  range  gate  in  the  time  intervcil, 
equals  — -^ • 

An  interfering  pulse  is  considered  to  ring  two  adjacent  range 
gates.  Therefore,  the  probability  of  false  alarm  in  one  range  gate  is 
not  stochastically  independent  of  the  time  of  occurrence  of  interfering 
pulses  in  either  adjacent  gate.  If  a  false  alarm  occurs  in  one 
range  gate,  the  probability  of  a  false  alarm  in  an  adjacent  gate  is 
increased.  The  effect  of  this  dependence  is  examined  in  Appendix  I. 
It  is  shown  that  an  assumption  of  independence  gives  results  which 
are  very  close,  in  the  cases  of  interest  here,  to  those  that  would 
be  obtained  by  very  lengthy  calculations  which  would  result  from 
the  exact  consideration  of  dependence.  Any  error  acts  to  make  the 
curves  presented  here  conservative  representations  of  the  exact 
results. 

Therefore,  T 

Pfa(T)  =  l-e   ^  =  1-6   ^]N   M 

For  false  alarm  time,  Tfa  >  1  ~  Pfa  ~  ^  (^^^ 


t  =  e 


i  -  /^N""^"^-f<^ 
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Tf^:  (In  2)  ('-— ^O  (23) 

fa 
now,  n where  n  is  the  false  alarm  number. 


n 


so,  P„  =  n^  ^  (2U) 


Equation  (23)  relates  pw  to  a  false  alarm  time,  while  equation  (2U) 
relates  p.,  to  a  false  alarm  number. 

The  equations  for  p  that  have  been  developed  so  far  have 
considered  only  the  effects  of  random  interfering  pulses.  To 
complete  the  picture,  the  effect  of  receiver  noise  must  be  considered. 
From  equation  (2),  the  probability  that  one  noise  pulse  exceeds  the 
established  bias  level  is 

Po  =  e  ^^  (25) 

Then  the  probability  that  an  interfering  pulse  or,  in  the  absence  of 
such  a  pulse,  a  noise  pulse  will  occur  in  a  chosen  range  gate  is     '-'' 

P(l)  =  ®   ^^  -  27j^S)-V-  2t{_S  (26) 

=  e  ^^(1  -  2$)^    2$ 

And  the  probability  of  a  false  alarm  in  the  chosen  range  gate  is 

N 


Figures  U2  through  53  are  further  working  curves  based  on  equations 
(26)  and  (2U).  With  a  given  density  of  interfering  pulses  and  a  given 
false  alarm  number,  the  resulting  bias  level  is  shown  for  various 
values  of  m  and  N» 

As  shown  in  the  equations  above,  the  bias  level,  y^,  is  determined 
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by  the  level  of  interference  and  the  desired  false  alarm  time.  Once 
this  bias  level  has  been  established,  it  is  possible  to  proceed  to  the 
calculation  of  the  probability  of  detection  of  a  target  by  the  radar. 
From  equation  (7),  the  probability  that  one  signal-plus-noise  pulse, 
with  signal-to-noise  ratio  X>  exceeds  the  bias  level  Y^,  is 

P^r  J  e'^^"^^^      lo  (2,rxf^)  dY  (2b) 

Figure  5U  shows  this  relationship  for  various  values  of  Y^' 

Using  this  single  signal  probability  p„,  the  desired  probability 
of  detection  is  the  partial  sum  of  the  binomial  expansion  from  M  to  N» 

H     „    i      N-i 

in  analyzing  the  double  threshold  receiver  the  false  alarm  number, 

7  K 

n,  was  chosen  as  6  x  10'   or  6  x  10  ,  corresponding  to  false  alarm 
times  of  one  minute  or  ten  minutes,  respectively,  for  a  radar  receiver 
having  a  one  microsecond  gate  width,  with  the  bias  level  thus 
established,  it  is  possible  to  display  the  relationship  between  the 
signal-to-noise  ratio,  X>  required  to  give  a  stated  probability 
of  detection  and  the  value  of  the  second  threshold  m»  pig^ires  2 
through  10  show  this  relationship  for  various  values  of  p^  and  N» 

pi  each  figure  there  is  evident  an  optimum  value  for  m»  This 
optimum  m  is  a  function  of  the  number  of  interfering  pulses  present; 
it  increases,  as  would  be  expected,  as  the  interference  density 
increases,  when  interference  is  present  the  curves  also  show  a 
minimum  value  for  u   below  which  the  desired  false  alana  time  cannot  be 
realized. 

In  each  calculation  all  of  the  interfering  pulses  were  assumed 
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to  be  large  enough  to  pass  the  first  threshold,  Yb»  The  increase 
in  signal-to-noise  ratio  required  for  a  constant  probability  of 
detection  as  the  numbers  of  interfering  pulses  increase  is  caused 
by  the  necessity  for  raising  the  first  threshold  to  reduce  the 
effect  of  the  receiver  noise.  The  probability  of  false  alarms  frcwi 
noise  must  be  reduced  in  this  case  in  order  to  maintain  the 
specified  over-all  false  alarm  time. 

on  each  curve  the  signal-to-noise  ratio  required  for  a  single 
threshold  integrating  receiver  in  the  absence  of  interfering  pulses 
is  shown.  It  should  be  remembered  that  the  specified  false  alarm 
time  applies  to  this  case  only  in  the  absence  of  interfering  pulses. 
It  is  not  possible  to  achieve  the  desired  T^^  with  this  type  of 
receiver  when  interfering  pulses  are  present  without  such  a  large 
increase  in  bias  level  that  practically  all  desired  echo  pulses  would 
be  eliminated. 

In  all  these  figures,  the  curves  relating  to  the  integrating  re- 
ceiver have  been  derived  from  conclusions  in  reference  (2). 

In  plotting  these  curves,  the  bias  level,  Yvj>  and  the  second  thres- 
hold, M>  were  chosen  so  as  to  achieve  the  given  false  alarm  time  in  the 
presence  of  interference.  The  signal-to-noise  ratio  required  for  the 
stated  detection  probability  was  calculated  in  the  absence  of  inter- 
ference. NO  attempt  was  made  to  consider  the  increased  probability 
of  detection  resulting  from  the  simultaneous  reception  of  an  interfering 
pulse  and  a  target  echo  pulse.  Therefore,  the  curves  represent  upper 
limits  for  X  in  the  presence  of  interference  and  give  a  conservative 
representation  of  the  results. 

The  expected  range  performance  of  a  radar  receiver  using  the 
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double  threshold  technique  may  be  seen  by  plotting  the  probability 
of  detection  in  one  scan,  p^^,  against  normalized  range,  ^/■dq* 
Rq  in  this  expression  is  the  solution  to  the  familiar  radar  range 
equation  for  a  signal-to-noise  ratio  of  one.  And 


R   _   1 


(30) 


Ro    (x)i 

Figures  11,  12,  and  13  show  these  range  curves  for  radar  receivers 
having  N  =  10,  N  =  20,  and  N  =  30.   In  each  figiire  range  performance 
cvirves  are  given  for  three  cases:  (1)  integrating  receiver  with  no 
interference,  (2)  double  threshold  receiver  with  no  interference,  and 
(3)  double  threshold  receiver  with  a  maximum  level  of  interference. 
The  double  threshold  receiver  in  conditions  of  no  interference 
exhibits  a  range  capability  approximately  i>^  below  that  for  an 
integrating  receiver  operating  under  the  same  conditions,  when 
interference  is  present  the  integrating  receiver  is  no  longer  able  to 
meet  the  false  alcirm  specifications.  The  double  threshold  receiver 
which  is  able  to  achieve  the  desired  false  alarm  time  has  a  range 
capability  between  $%   and  10^  below  that  for  the  integrating  receiver 
operating  with  no  interference. 

The  range  capability  of  the  double  threshold  receiver  is 
relatively  insensitive  to  changes  in  the  desired  false  alarm  time. 
Figures  lU  and  15  compare  the  range  performances  for  several  valiies 
of  false  alarm  number,  n.  For  a  probability  of  detection  of  one-half 
an  increase  in  false  alarm  number  from  10^  to  10^^  (corresponding  to 
an  increase  in  false  alarm  time  from  0.01  seconds  to  nearly  3  hours 
for  a  one  microsecond  range  gate  radar  receiver)  results  in  a  loss 
of  maximum  range  capability  of  about  20^. 
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CHAPTER  IV 
THE  DOUBLE  THRESHOLD  RFCEIVER  -  ^FLUCTUATING  TARGET  ECHOES 

The  expected  performance  of  the  double  threshold  radar  receiver 
with  target  echo  signals  which  are  fluctuating  in  amplitude  can  be 
calculated  with  the  results  of  the  previous  section  as  a  basis. 
Rather  detailed  pulse  studies  of  the  target  echo  signals  from  air- 
craft indicate  that  the  most  reasonable  probability  density  for 
the  signal-to-noise  power  ratio  is 

(X,  x)  ^ 

where  X  =  signal-to-noise  power  ratio 

Y  •  average  of  X  over  all  fluctuations . 
This  density  must  be  applied  to  the  case  where  the  target  echo 
amplitude  fluctuates  independently  from  pulse  to  pulse,  and  must 
be  applied  to  the  case  where  the  target  amplitude  is  assumed  to 
be  constant  for  a  single  scan  (n  periods)  but  to  fluctuate  indepen- 
dently from  scan  to  scan.  Most  actual  targets  would  probably 
present  fluctuations  which  were  somewhere  between  these  two  extreme 
situations. 

From  equations  (9)  and  (11),  the  probability  that  one  signal-plus- 
noise  pulse  with  average  signal-to-noise  power  ratio  x  exceeds 
the  bias  level  y^j  is 

por  pulse  to  pulse  fluctuations,  the  probability  of  detection  in  a 
single  scan  with  a  double  threshold  receiver  is 
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N 

i=M   ^ 
As  with  the  case  of  the  non-fluctuating  target  echoes,  the  bias  level, 
Yvj,  is  determined  by  the  desired  false  alarm  time.  This  time  is  not 
affected  by  any  fluctuations  of  the  target  echo  signal. 

Figures  16  through  21  show  the  signal-to-noise  ratio  as  a  function 
of  the  second  threshold  when  pulse  to  pulse  fluctuations  are  assumed. 
For  canparison  purposes  the  same  gate  width,  false  alarm  time,  and 
numbers  of  interfering  pulses  have  been  chosen  as  for  the  curves  in  the 
corresponding  non-fluctuating  cases.  The  curves  have  the  same  general 
shape  as  for  the  curves  in  the  non-fluctuating  situation.  Again  there 
is  a  definite  optimxan  value  for  m  and,  when  interference  is  present, 
there  is  a  definite  minimum  value  for  u.  This  minimum  value  of  the 
second  threshold  is  the  same  as  for  the  corresponding  non-fluctuating 
cases.  The  pulse  to  pulse  target  fluctuations  cause  an  increase  in  the 
signal-to-noise  ratio  that  is  required  to  achieve  a  stated  probability 
of  detection  for  values  of  u     approaching  N*  as  compared  to  the  non- 
fluctuating  target  echoes.  Note  that  this  increase  made  it  necessary 
to  use  a  two  cycle  logarithmic  scale  for  X,  while  a  one  cycle  logarith- 
mic scale  has  been  adequate  for  all  previous  figures. 

Figures  22  through  2U  show  the  probability  of  detection  as  a 
function  of  the  normalized  range  yihen   pulse  to  pulse  target  fluctuations 
are  considered.  These  curves  correspond  to  figures  11  through  13  for 
the  non-fluctuating  case.  A  canparison  of  the  three  figures  indicates 
that  the  optimum  range  performance  relative  to  an  integrating  radar 
receiver  requires  use  of  high  values  of  n  when  pulse  to  pulse  fluctu- 
ations are  present.  For  N  =  30  the  double  threshold  receiver  with  no 
pulse  interference  has  a  maximum  range  loss  of  about  lOjg  relative  to  an 
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integrating  receiver.  For  high  interference  levels  this  double 
threshold  receiver  has  a  range  loss  of  about  20^^  relative  to  an 
integrating  receiver  operating  without  interference. 

When  scan  to  scan  target  fluctuations  are  considered,  the  ex- 
pression for  the  probability  of  detection  becomes 

(^  -x/x- 
Pd=0-rr-e    Pd'  dx  (3U) 

0  X 
where  p^j'  is  the  detection  probability  for  the  corresponding  non- 
fluctuating  case.  This  expression  cannot  be  solved  in  closed 
form  and  graphical  integration  methods  have  been  used  to  obtain  p^. 
Figures  25  through  27  show  the  probability  of  detection  as  a  function 
of  the  normalized  range  for  scan  to  scan  fluctuations.  The  range 
performance  of  the  double  threshold  radar  receiver  relative  to  the 
integrating  receiver  is  seen  to  correspond  approximately  to  the  non- 
fluctuating  case.  A  range  loss  of  $%   to  10^  occurs  in  the  double 
threshold  receiver  as  compared  to  the  interference-free  integrating 
receiver.  This  range  performance  is  seen  to  be  apparently  independ- 
ent of  N. 

A  comparison  of  the  fluctuating  and  non-fluctuating  target  echo 
cases  is  shown  in  Figure  28.  For  reasonably  high  probabilities 
of  detection,  target  fluctuation  is  seen  to  have  a  very  great  effect 
on  range  performance.  These  cujrves  correspond  in  shape  and  in  rela- 
tive position  to  those  computed  by  sterling  for  the  integrating 
radar  receiver 


a 
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CHAPTER  V 
ANTENNA  BEAM  SHAPE  EFFECTS 

The  effect  of  the  beam  shape  as  a  radar  antenna  scans  past  the 
target  can  be  considered  as  a  pulse  to  pulse  fluctuation  with  a 
probability  density  function  determined  by  the  squared  power 
pattern  (the  two-way  pattern).  A  reasonable  assumption  for  the 
beam  shape  is  » 

f  (©)  =  X  =  Xin  e         for  X^  iXa     (35) 
where  3Qa  is  the  signal-to-noise  power  ratio  at  the  center  of  the 
beam,  l^e  details  of  the  derivation  of  the  beam  density  function  are 
given  in  Appendix  li.  The  resulting  beam  density  function  is 

O.U23  ^    ^ 

Pb  = pr::zr3':-.    ^Xm  ^  x  et  Xo     (36) 

X  ^jln  xa/x^ 


Graphical  integration  methods  were  used  to  determine  the  probability 
of  detection.  Figures  29  through  32  show  the  x  vs.  M  curves  for 
N  ■  10  and  N»  20,  comparing  the  curves  obtained  in  the  non-fluctuating 
case  with  those  obtained  by  considering  beam  shape  losses. 

Consideration  of  beam  shape  does  not  alter  the  general  shape 
of  the  curves.  It  introduces  a  constant  loss  in  all  curves.  The 
effect  of  beam  shape  can  be  included  in  range  calculations  by 
applying  a  1.5  db  correction  to  the  signal-to-noise  power  ratio. 
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CHAPTER  VI 
CONCLUSIONS 

Analysis  of  the  figures  of  performance  of  the  double  threshold 
type  of  radar  receiver  lead  to  the  following  general  conclusions: 

1.  The  double  threshold  technique  permits  the  retention  of 
any  reasonable  false  alarm  time  in  the  presence  of  large  numbers 
of  interfering  pulses. 

2.  Reasonable  false  alarm  times  are  not  possible  ?rith  an  inte- 
grating radar  receiver  in  the  presence  of  large  numbers  of  inter- 
fering pulses.  The  double  threshold  receiver  which  does  give  the 
desired  false  alarm  time  results  in  a  loss  of  maximiam  radar  range 
of  from  $%   to  10^  as  compared  to  the  integrating  receiver  operating 
in  the  absence  of  interference. 

3.  TJ-or  any  value  of  n>  the  number  of  target  echoes  vrfiich  are 
received  when  an  antenna  scans  past  the  target,  there  is  an  optimum 
value  of  M>  the  second  threshold  number.  This  optimum  M  increases 
with  increase  in  the  number  of  interfering  pulses  which  must  be  handled. 

U*  For  any  number  of  interfering  pulses  there  is  a  minimum  value 
of  M  below  which  the  desired  false  alarm  time  cannot  be  realized. 

5.  Fluctuation  in  the  amplitude  of  target  echoes  results  in  minor 
modifications  in  the  value  of  M  which  is  optimum,  but  fluctuation  does 
not  alter  the  minimum  value  of  M  which  permits  desired  operation. 
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APPENDIX  I 
STOCHASTIC  DEPENDENCE  OF  FALSE  AURMS  IN  ADJACENT  RANGE  GATES 


The  probability  of  at  least  one  false  alarm  in  time  T  is : 


Pfa(T)  =  1  - 


(l  -  Pmi)  (  ^  -  Pm2)(^  -  Pm3) (^  - 


where  p.,  is  the  probability  that  pulses  will  occur  in  the  n-th  range 
gate  on  M  or  more  of  N  successive  repetition  periods  (a  false  alarm 
in  the  n-th  range  gate);  and  r,  the  last  range  gate  in  the  time 
interval  T,   equals   "   . 

The  probability  of  false  alarm  in  one  range  gate  is  not  independent 
of  the  time  of  occurrence  of  interfering  pulses  in  either  adjacent 
range  gate.   If  a  false  alarm  occurs  in  one  range  gate,  the  probability 
of  a  false  alarm  in  an  adjacent  gate  is  increased.  However,  as  only 
alternate  gates  are  considered,  the  probabilities  of  false  alarms 
are  stochastically  independent. 

for  even  numbered  gates  -  -  P^^^^^^  -  1  -  ^^l-p^^)  (l-p^^ J  -  -  -J 

r 
=  1  -  e 

r  ^      "1 

for  odd  numbered  gates P-     =  1  -  (l-Pw  )(l~Pf,  ) 1 


f  % 


=  1  -  e 
''fa  C^)  =  P(fa     or  fa  ^ 


—    p  \         p        _    p 

^^even"^  ^^odd    ( ^^g^en  ^^  ^^odd )  • 
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(fa_en  and  fa  ,,)  "  ^^even      ( ^^o^i^i  f agven ) 


so,  Pfa  (T)  =  2  (1-e  2^")  -  (1-e   ")  (1-e  ^  ^^^  ) 

v»here  (/i    is  the  factor  showing  the  increased  probability  of  false 

alarm  in  one  range  gate  when  there  has  been  a  false  alarm  in  an 

adjacent  gate. 

-rp  -pp 

P^_  (T)  =  1  -  e   M  _^  (1  _  e  "^  ^) 


r 


■2Pm 


fa    ^  =0-  -  ®    +^'>  ) 
even  ^      ^ 

then  0  <<A  -C   ®    • 
For  the  values  of  m»  N>   and  jx.     are  encoxintered  in  this  discussion 
of  the  double  threshold  receiver,^  was  found  to  be  ver>  close  to 
zero.  Typical  values,  obtained  by  enumerating  the  ways  of  obtaining 
a  false  alarm  in  one  range  gate  when  a  false  alarm  has  occurred  in  the 
adjacent  gate,  are  : 

10  5  0.35 

10  7  0.077 

20  10  0.115 

20  lU  0.00U3 

By  assuming  ^  =  0,  the  curves  of  false  alarm  number  and  of  the  per- 
formance of  the  double  threshold  receiver  are  conservative  represen- 
tations of  the  exact  values,   in  most  cases  the  difference  between 
the  curves  obtained  by  assuming  c^    z  0   and  by  calculating  the 
correct  value  of  (/^     is  too  small  to  be  shown  graphically. 


2ti 


Therefore , 

=  1  -  6-"-%  :  1  -  e"  ^A^P«. 


Pf^(T) 


when  T  equals  the  false  alarm  time,  Tfg.  >  then 
1  -  Pf  a  =  i 


and       -^  =  e  '^ 
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APPENDIX  II 
DERIVATION  OF  ANTENNA  BEAM  SHAPE  DENSITY  FUNCTION 


Q   .,.-> 


AG^' 


Ax; 

4  Iq) 


.  AS 


AX 


It  is  reasonable  to  assume    ^  ( 0  )  =  X  = 


^< 


or 


Y  ^  M 


-cX©^ 


^ 


^°^  X  ^  J  Xm. 
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since  the  same  antenna  is  assumed  for  transmitted  and  received  signals, 
the  bearawidth  is  calculated  between  the  ^X^  points  of  this  two-way 
power  beam  pattern. 

■f  (©)  "  -2.xeXMe  -   -2=K©x 


-f'(9)--   -2<^.X-  J^-  In  J^ 


a 


r.i*.  o»»Mr '(*»^  «3 


X 


7- 


e  =  i.ig 


X 

The  probability  that  a  single  signal-plus-noise  pulse  will  exceed 
the  bias  level  Y^  is      f^^ 


In  this  equation,  p  is  the  probability  of  detection  in  one  range  gate 

A 

if  beam  shape  is  not  considered.  The  solution  of  this  beam  shape 
equation  for  py«  has  been  carried  out  by  graphical  integration. 
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